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Hot dust in normal star-forming galaxies: JHKL' photometry 
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Abstract. We present JHK and 3.8 um (L r ) photometry of 26 galaxies in the Infrared Space Observatory (ISO) 
Normal Galaxy Key Project (KP) sample and of seven normal ellipticals with the aim of investigating the origin 
of the 4/im emission. The majority of the KP galaxies, and all the ellipticals, have K — L £ 1.0, consistent with 
stellar photospheres plus moderate dust extinction. Ten of the 26 KP galaxies have K — L 1.0, corresponding 
to a flat or rising 4 /xm continuum, consistent with significant emission from hot dust at 600-1000 K. K — L is 
anticorrelated with ISO flux ratio -F6.75/-F15, weakly correlated with line ratio [Oi]/[Cii], but not with [Cn]/FIR 
or IRAS ratio Jeo/fioo. Photodissociation-region models for these galaxies show that the hot dust responsible 
for red K — L resides in regions of high pressure and intense far-ultraviolet radiation field. Taken together, these 
results suggest that star formation in normal star- forming galaxies can assume two basic forms: an "active", 
relatively rare, mode characterized by hot dust, suppressed Aromatic Features in Emission (AFEs), high pressure, 
and intense radiation field; and the more common "passive" mode that occurs under more quiescent physical 
conditions, with AFEs, and without hot dust. The occurrence of these modes appears to only weakly depend on 
the star-formation rate per unit area. Passive star formation over large scales makes up the bulk of star-forming 
activity locally, while the "active" regime may dominate at high redshifts. 

Key words. Galaxies: spiral - Galaxies: starbursts - Galaxies: ISM - ISM: Dust 



1. Introduction 

The Infrared Space Observatory (ISO) mission has pro- 
vided an unprecedented view of the interstellar medium 
(ISM) in galaxies from the near-infrared (NIR) contin- 
uum between 3 and 5 /im to the C + fine-structure transi- 
tion at 158 /im and beyond. Infrared spectra (3 to 12 /im) 
obtained by ISO-PHOT reveal the mid-infrared (MIR) 
emission from normal galaxies to be characterized by the 
Aromatic Features in Emission (AFEs) at 3.3, 6.2, 7.7, 
8.6, and 11.3 /im, and an underlying continuum which con- 
tributes about half the luminosity in the 3 to 12 /mi range 



Hydrogenated Carbon particles (HACs), and their rela- 
tive contribution in different environments is still under- 



way (Jenniskens & Desert 1993; Lu 1998; Cesarsky et al. 
|1998[ ). Mor eover the relationship between the carriers of 
the AFEs and the underlying continuum remains obscure 
(Boulanger et al. |l996j Helou et al. [2000D . 



(Helou et al. |2000p . 

However, the ISO data have also raised many 
questions about the ISM constituents. In particular, 
the nature of the dust responsible for the AFEs is 
not yet clear, and the debate over three-dimensional 
Very Small Grains (VSGs), two-dimensional Polycyclic 
Aromatic Hydrocarbon molecules (PAHs), or Amorphous 



A related question is the relative contributions of dust 
in the ISM, stellar photospheres, and circumstellar dust 
to the MIR radiation in spirals and ellipticals. IRAS data 
already convincingly showed that the MIR emission in spi- 
rals is mainly due to a continuum from small grains tran- 
siently heated to high temperatures and AFEs. On the 
other hand, in elliptical galaxies the ISM-to-stellar ratio 
is generally low, so that photospheric and circumstellar 
emission from evolved red giant s may contribute signifi- 
cantly to the MIR (Knapp et al. |1992| , Mazzei & de Zotti 



1994 ). It has been argued, tho ugh, that the same IRAS 
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color-color relation (e.g., Helou 198E) holds for ellipticals 
and spirals, which implies a similar ISM origin (Sauvage 
& Thuan 



1994D 



TIRGO, Gornergrat, 



The continuum component between 3 and 5 /im can 
be extrapolated fairly well to the continuum level at 9- 
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10 /im (Helou et al. 2000 ), which suggests that some of 
the 3-5 /tm continuum must be produced by the ISM. 
Nevertheless, a portion of the 3-5 /im continuum must 
arise from stellar photospheres since they dominate the 
emission at 1-2 /im. Although the signal-to-noise in the 
3-5 /jm ISO spectra is relatively low, that is the spectral 
region where the MIR spectra of the 43 galaxies observed 
by Helou et al. ( |2000 ) show the most significant galaxy-to- 
galaxy variation. To understand the physical origin of this, 
and assess the relationship of the MIR continuum to the 
ubiquitous AFEs, it is necessary to separate the ISM con- 
tribution from that of stellar photospheres and circumstel- 
lar dust. To this end, we have acquired JH KL' p hotom- 
etry of the galaxies studied in Helou et al. ( 2000| ), which 
constitutes one of the ISO Normal Galaxy Key Projects. 
A variety of ISO observations have been secured for this 
sample, including mid-infrared spectra (Helou et al. 2000] ), 
ionic and atomic fine-structure line fluxes (Malhotra et al. 
1997, Malhotra et al. 2001), mid- infrared maps at 7 and 



15 /im, and a small number of 4.5 /im images (Dale et al. 
2000| ). 

The rest of the paper is organized as follows: Sect. 
H describes the observations, reduction, and photometric 
calibration. We compare the JHKL colors of the sample 
galaxies with their ISO properties in Sect. ||, and with 
photodissociation-region models in Sect. ||. Section || dis- 
cusses the nature of the 4 /tm continuum, and our interpre- 
tation of these results in terms of "active" and "passive" 
star formation. 

2. The sample, the observations, and the colors 

The ISO Key Project for Normal Galaxies (PI: G. Helou, 
KP) was proposed to study the ISM of a broad range 
of normal galaxies using several instruments aboa rd ISO. 
The sample, described in Malhotra et al. ( 2001 ), is de- 
signed to capture the great diversity among galaxies, in 
terms of morphology, luminosity, far-infrared-to-blue lu- 
minosity ratio FIR/B, and IRAS colors. The project ob- 
tained ISO observations of 69 galaxies, including nine rel- 
atively nearby and extended objects. The remaining 60 
galaxies cover the full range of observed morphologies, lu- 
minosities, and star-formation rates seen in normal galax- 
ies, and includes five dwarf irregulars discussed in Hunter 
et al. ( |2001| ). Two of these, NGC 1156 and NGC 1569, are 
part of our observed sample, and are also the two nearest 
objects (6.4 and 0.9 Mpc, respectively). The median dis- 
tance of the KP sample (including these nearby objects, 
and with H = 75 km/s/Mpc) is 34 Mpc (Dale et al. |2000D . 
This means that with the 14" observing aperture we are 
sampling a central region roughly 1 kpc in radius. We ob- 
served 26 KP galaxies. 

We also selected seven quiescent ellipticals from the 
sample of early-type galaxies observed with ISO by 



2.1. Observations 

The observations were acquired at the 1.5 m f/20 Infrared 
Telescope at Gornergrat (TIRGOQ), with a single-element 
InSb detector. The photometer is equipped with stan- 
dard broadband filters (J 1.2 /im, H 1.6 /tm, K 2.2 /tm, 
and L' 3.8 /tm) with diaphragms in the focal plane defin- 
ing the aperture dimensions. The galaxy coordinates were 
taken from NEDQ and checked with the DSS^]; all observa- 
tions were acquired with a 14" aperture, after maximizing 
the infrared signal. Sky subtraction was performed with 
a wobbling secondary at frequencies that ranged from 2.1 
to 12.5 Hz, according to the integration time for the in- 
dividual measurement. The modulation direction was in 
a EW direction, with an amplitude of roughly 3'. Beam 
switching was used to eliminate linear variations in sky 
emission. 

Photometric calibration was achieved by similarly ob- 
serving several standard stars nightly from the CIT (Elias 
et al. [1982D and the ARNICA (Hunt et al. |l998j ) standard 
lists. Nightly scatter of the photometric zero point was 
typically 5% or better in JHK and 8% or better in L' . 

The final photometry for the 33 galaxies observed is re- 
ported, together with their basic data, in Table [|. In what 
follows, we have transformed K—L' to K — L (this last is 
center ed at 3.5 /tm) using the transformation by Bessell & 
Brett ( [1988 ) , and we will use K — L to denote such colors. 
In the various plots and when testing for correlations, the 
NIR data have also been correct ed for i) Galactic extinc- 
t ion a ccording to Schlegel et al. ( [1998 ) and Cardelli et al. 
( 1989| ), and ii) K dimming using the precepts mentioned 
in Hunt & Giovanardi ( |1992| ). 



2.2. JHKL colors 

The NIR colors of the observed galaxies are shown in 
Fig. 0. The filled symbols represent those (10) galaxies 
with K - L £ 1.0, that is a rising continuum (see be- 
low). Also shown in the diagra m are Hn galaxies taken 
from Glass & Moorw ood ( 1985 ) and starbursts from Hunt 
& Giovanardi ( |l992 ). Three low-metallicity blue compact 
dwarfs (BCDs) are also s hown as asterisks: NGC 5253 
(Glass & Moorwood |l985| ), IIZw40 (Thuan |l983|) , but 
no J), and SBS 0335-052 (Hunt et al. |200l| ). The BCDs 
have the most extreme colors since they tend to be quite 
blue in J — H (because of low metallicity and youth), and 
red in K — L (because of ionized gas and hot dust). 

Mixing curves show how the colors change when vari- 
ous physical processes increasingly contribute to the emis- 
sion observed (see Hunt & Giovanardi 1992] ); they illus- 



Malhotra et al. (200C), or with published JHK photome- 
try (Frogel et al. 1978). Several of the latter do not have 
the rich array of ISO observations obtained for the KP 
galaxies. 



1 The Infrared Telescope at Gornergrat (Switzerland) is op- 
erated by CAISMI-CNR, Arcetri, Firenze. 

2 The NASA/IPAC Extragalactic Database (NED) is oper- 
ated by the Jet Propulsion Laboratory, California Institute of 
Technology, under contract with the U.S. National Aeronautics 
and Space Administration. 

3 The DSS was produced at the Space Telescope Science 
Institute under U.S. Government grant NAG W-2166. 



Table 1. Basic data and NIR photometry 
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Basic data are from NED. The errors reported for the NIR photometry (columns 8 through 11) are la. 
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Fig. 1. NIR color-color diagram; the K — L color in the 
top panel is transformed from K — L' as described in the 
text. Filled symbols refer to galaxies with K — L > 1.0. 
The 7 quiescent ellipticals are indicated with an E. Shown 
as sta rs are Hn galaxies taken from Glass & Moor wood 
( 1985 ) and starbursts from Hunt & Giovanardi ( 1992 ), and 
shown as asterisks the BCDs (see text). Mixing curves 
show how the colors change when various physical pro- 
cesses increasingly contribute to the emission observed 
(see Hunt & Giovanardi 1992); the end points of the mix- 



ing curves indicate equal if-band contributions from stars 
and the process (hot dust, ionized gas, A stars). The tick 
marks on the extinction line show unit Ay magnitude in- 
crements. 



trate that, on the basis of NIR color, it is possible to dis- 
tinguish among stellar photospheric emission, "passive" 
dust in extinction, and "active" dust in emission. H — K 
tends to be red (> 0.35) for both dust extinction and 
dust emission, while K — L is red only because of emis- 
sion by dust or ionized gas (see the free- free line). Thus 
red K — L can be used to signal a substantial contribution 
from hot dust emission in the observing aperture. The 600 
and 1000 K mixing curves span the observed colors quite 
well, and although the VSG emission is not thermalized, 
the curves show that the NIR colors are consistent with 
dust at these temperatures. 

The K — L color can also provide an estimate of 
the slope of the spectral energy distribution (SED) be- 
tween 2 and 4/im. Using the zero-magnitude fluxes given 



Therefore^, K — L 1.0 is where the slope a changes 
sign and becomes negative, signifying a continuum /„ ris- 
ing with increasing wavelength; when K — L < 1.0, a is 
positive. The mean a averaged over the entire KP sample 
is > 0, with /„ oc v +0M5 (Helou et al. |2000p. If we include 



the ellipticals, the median K — Loi our 33 galaxies is 0.50, 
which gives a — 0.90, steeper than, but consistent with, 



Helou et al. (gOOOJ). If only the ISO Key Project galax- 
ies are considered, median K — L = 0.68, corresponding 
to a = 0.54, remarkably close to that reported in Helou 
et al. The median K — L for the ellipticals only is 0.21, 
which corresponds to a very steep falling continuum with 
a = 1.5. 

The mixing curves shown in Fig. [l] illustrate what frac- 
tion of the observed flux is due to hot dust; they assume 
a mix of stellar photospheres with intrinsic stellar color 
(K — L)» = 0.3, or (K — £/)„ = 0.5, plus dust emis- 
sion. When the SED is flat at 4 [im(K — L w 1.0, slope 
a = 0), 600 K hot dust comprises roughly 5% of the total 
.ftT-band flux; hotter dust (e.g., 1000 K) would constitute 
30%. It is unlikely that flat or rising continua are due to 
free-free emission from ionized gas, since even with a 50% 
emission fraction from gas, the continuum is still falling 
\(K — L) = 0.7]. This means that K — L, or alternatively 
the slope of the 4 /xm continuum, is a remarkably sensitive 
diagnostic of hot dust: small fractions (5-30%) of dust 
emission cause large variations (0.6 mag) in the K — L 
color. A 50/50 if-band mix of hot dust and stars would 
produce a K — L color between 1.2 and 2.7, depending 
on the dust temperature; NGC4519, the galaxy with the 
reddest H — K (= 0.82), may contain such a mix. 

3. Comparison with ISO photometry and spectra 

If the non-stellar 3.8 /im flux is truly associated with hot 
dust, we would expect to find correlations between ISO 
MIR observations and L', if both arise from a similar 
small-grain population. Also, observations of the two most 
important photo-dissociation region (PDR) cooling lines, 
[C n] and [O i] , suggest that gas and dust temperatures in- 
crease together (Malhotra et al. 2001); we might therefore 
expect trends with these line fluxes and ratios and K — L. 

3.1. Broadband fluxes 

We first checked that our ground-based (small-aperture) 
L' fluxes were broadly consistent with the ISO fluxes at 
-^6.75 an d Fi5 (Dale et al. 200C). All three are mutually 
correlated (not shown), except for two cases (NGC278 
and NGC 5713) where the ISO flux is substantially larger 
than expected from the general trend. This discrepancy is 
almost certainly due to aperture effects, since the ground- 
based data are acquired in a 14" aperture, while the ISO 
values published in Dale et al. ( 2000] ) are total fluxes as 
extrapolated from curves of growth. NGC 1569 is gener- 



in Koornneef (1983), and assuming a power-law depen- 
dence f v oc v a , we find that K — L = —0.50 a + 0.95. 



4 The intercept climbs from .95 to 1.10 if we adopt instead 
the scale of Wilson et al. ((L972|) . 
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Fig. 2. Hybrid ground-based/ISO color-color plot: the top 
panel shows the ratio of the 3.8 /im flux (not transformed 
to L) and F±s plotted against K — L, and the bottom panel 
the ratio of ^6.75 and F±5 plotted against K — L. The 
dotted line shows the best-fit regression; the correlation 
coefficients are shown in each panel (parentheses denote 
the coefficient without the two nearest galaxies NGC 1569 
and NGC 1156). Solid dots are used when the corrected 
H — K > 0.35. 



ally an outlier; however, its properties are similar to those 
of other low-metallicity dwarf irregulars (see Hunter et al. 
2001). The other dwarf irregular in our observed sample is 
NGC 1156, again a clear outlier in most of the subsequent 
plots and correlations. As noted in Sec. ||, these galaxies 
are also the two closest galaxies, which considerably com- 
plicates the comparison of the ISO and our ground-based 
14" aperture data. 

We next compare ISO colors with K — L in Fig. || 
where F 6 ,r 5 /Fi 5 and the "hybrid" color i^.s/^is are plot- 
ted against K—L. Ground-based K—L turns out to be cor- 
related (2.5 a) with Fq, 75/^15^, but is uncorrelated with 
the hybrid color. The sense of the ^6.75/^15 vs. K — L cor- 
relation is such that redder K — L implies lower F§ ^/ F\§. 
NGC 1569 stands out since it has a very low F, g. 75/-F15 
ratio for its K — L color. As an extreme case, in Fig. 0, 
we have also plotted SBS 0335-052 (not considered in the 
correlation), using the ISO data from Thuan et al. (1999) 
and K-L from Hunt et al. (pOOT|). SBS 0335-052 is an 



unusual BCD with 1/40 solar metallicity, and the ISO 
spectrum shows no AFEs; not surprisingly Fig. ^ shows 

5 If the nearest objects, NGC 1569 and NGC 1156, are elim- 
inated from the regression the significance becomes 4.0 a. If, 
in addition, the reddest K — L value (IRAS 23365+3604) is 
eliminated, the significance becomes 2.1 a. 
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Fig. 3. IRAS 12fim/25fj,m flux ratio vs. K - L. No IRAS 
upper limits are shown. Solid points refer to galaxies with 
a corrected H — K > 0.35. The best-fit regression is shown 
together with the correlation coefficient (in parentheses is 
that without NGC 1156 and NGC 1569, the two nearest 
galaxies). 



that this galaxy has a lower F e 75/F15 ratio than any of 
the galaxies in the ISO sample. 

Since dust temperatures in large "classical" grains are 
connected with the IRAS flux ratio Fqq/Fiqq, and since 
this ratio and -Fe.75 /-P15 are anticorrelated (Dale et al. 
|2000D , we might expect K — L to also be anticorrelated 
with Fqq/Fiqq. We found no such correlation but note, 
instead, that IRAS -F12/F25 is anticorrelated (2.6 a) with 
K — as shown in Fig. Together with the trend with 
K — L and -F6.75/-F15, this means that the presence of 
hot dust is: a) usually linked to the suppression of AFEs 
which dominate the 12 ^m band (Helou et al. |1991| ), and 
b ) largely independent of the temperature and character- 
istics of the large grains. 



3.2. Line measurements 

In PDRs and H11 regions, gas and dust are intimately 
related. In ionized regions, gas and dust compete for far- 
ultraviolet (FUV) photons, but the neutral gas in PDRs is 
heated predominantly by photoelcctrons from small dust 
grains (Watson [l972| ; Hollenbach & Tielens |l997| ) . Neutral 



6 If the object with the reddest K — L (IRAS 23365+3604) is 
eliminated from the regression, the significance becomes 2.4 a. 

7 Although these galaxies were selected to be sufficiently dis- 
tant, given the size of the IRAS beams compared to our ground- 
based photometry, aperture effects almost certainly contribute 
significantly to the scatter. 
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Fig. 4. ISO emission-line ratios vs. K — L. The top 
panel shows the log of the ratio of [C n] (158 /zm) and 
[Oi] (63 fim) line fluxes vs. K — L, and the bottom panel 
the log of the ratio of [C n] and FIR vs. K — L. The corre- 
lation coefficients are shown in each panel (the correlation 
coefficient in parentheses is that without NGC1156 and 
NGC 1569, the two nearest objects). Solid dots refer to 
galaxies with corrected H — K > 0.35. 



gas is cooled primarily by atomic and ionic fine-structure 
lines, [Cn] (158 /im) and [Oi] (63 /im), and these lines 
can be used as diagnostics for the physical conditions in 
the PDR gas (Tielens & Hollenbach |1985|) : the [Oi] line 
is expected to become more important relative to [C n] 
in warmer and denser gas. [Cn] and [Oi] line fluxes for 
the galaxies in our sample have been measured by ISO 
(Malhotra et al. 2001), and in this section we analyze those 
measurements in the context of our new photometry. 

An important diagnostic is the ratio of [C n] and the 
far- infrared flux (FIR), since it measures essentially the 
efficiency of the photoelectric heating of the gas by dust 
grain ejection. This ratio [C ii]/FIR tends to decrease with 
warmer FIR colors Fqq /Fioo and increasing star- formation 
activity as indicated by FIR/B (Malhotra et al. 2001 ). 
Moreover, warmer gas, as signified by smaller [Cii]/[0i], 
correlates with warm dust or larger Feo/Fioo; this last is 
the most significant correlation in the study by Malhotra 
et al 



(2001). Since red K — L should be related to hot 
dust and its temperature, we might expect to find similar 
correlations with normalized FIR line fluxes and ratios. 
However, this supposition is not borne out by the data 
(see Fig. ^): we find no correlation between K — L and 
[C n] /FIR, and only a weak anticorrelation between K — L 
and [Cii]/[Oi]. Again, K — L appears to be measuring a 
hot-dust phase, not closely connected to the properties of 
the cooler dust. 
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Fig. 5. Logn (upper panel) and Log Go (lower) plotted 
against K — L. The PDR parameters are taken from 



Malhotra et al. (2001). The dotted line in the lower panel 
shows the best-fit regression; the correlation coefficients 
are shown in each panel. Solid points refer to galaxies with 
corrected H — K > 0.35. Estimate of correlation coeffi- 
cients and regression include all objects (the parentheses 
contain the coefficient without NGC 1156 and NGC 1569, 
the two nearest galaxies). 

4. PDR models 

Photodissociation models for all the galaxie s in o ur sam- 
ple have been reported in Malhotra et al. ( 2001 ), calcu- 
lated from the model grids of Kaufman et al. ( 1999 ). These 
models calculate the line emission in [C n] and in [O i] (in- 
cluding the line at 145 /im) and the dust continuum emis- 
sion for a plane slab of gas illuminated on one side by FUV 
radiation. Gas heating in the models is dominated by pho- 
toelectrons ejected from classical (large) grains an d from 
VSGs according to the recipe of Bakes & Tielens ( 1994 ) . 
From these models, it is possible to infer values for the 
neutral gas density n and the FUV radiation field Go, as 
well as for the pressure P and temperature T of the gas. 

Figure || shows gas density n and FUV flux Go, as de- 
rived from PDR models, plotted against K — L. Without 
NGC 1156 (one of the two nearest galaxies), K — L is well 
correlated with Go, at the 3.5 a level. This is a clear indica- 
tion that red K — L is connected to strong FUV radiation 
fields. Also shown in Fig. | is SBS 0335-052, with a G 
~ 10000 times stronger than the local value (Dale et al. 
[2001 ); this galaxy follows the same correlation. 

We then examine whether the primary dependency of 
red K — L color is on neutral gas temperature T or on 
its pressure P, given that the K — L color is independent, 
or roughly so, of neutral gas density n. The relevant PDR 
parameters are plotted in Fig. where T and P as in- 
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Fig. 6. LogT (upper panel) and LogP (lower) plotted 
against K — L. The PDR parameters are taken from 
Malhotra et al. (2001). The dotted lines show the best- 
fit regression; the correlation coefficients are reported in 
each panel. Solid points refer to galaxies with corrected 
H — K > 0.35. Estimates of regression line and correlation 
coefficients include all galaxies (the parentheses contain 
the coefficient without the two nearest objects, NGC 1156 
and NGC 1569). 



ferred from the PDR models are plotted versus K — L. 
The thermal pressure P of the PDR is expected to be 
approximately equal to the thermal pressure of the adja- 
cent Hn region (Malhotra et al. 2001 ). The correlation 
between P and K — L is weak but becomes significant 



(2.8a) once we exclude NGC 1156. There is no apparent 
correlation between T and K — L. Taken together, these 
results suggest that the PDR pressure (roughly equal to 
the H II pressure) and the FUV radiation field Go are the 
main factors which govern the existence of hot dust in star- 
forming galaxies; neutral gas density and temperature by 
themselves are less important influences. 



5. The nature of the 4 /im continuum 

We first discuss the nature of the 4 /jm continuum, as in- 
ferred from our observations. Because of the filter cut- 
off, L' misses the 3.3 /im aromatic emission feature (e.g., 
Moorwood & Salinari 1983; Moorwood 1986), and thus 
should be a measure of the AFE-free continuum. As men- 
tioned in Sec. |2.2| , the K — L color is a sensitive diagnos- 
tic of the fraction of hot dust contained in the observing 
aperture. The 4 /im continuum of most of the galaxies we 
observed is consistent with stellar photospheres together 
with moderate dust extinction. As expected, the best cases 
for pure stellar K — L colors are the ellipticals which have 



median K — L = 0.21. Within the 14" observing aperture, 
we find no evidence for circumstellar dust in the ellipticals, 
since their NIR colors are inconsistent with dust emission 
at any temperature. 

In roughly a third of our sample, the K—L color ( }t 1) 
signifies a flat or rising 4 /im continuum. Even a small 
fraction (5%) of hot (600 K) dust can cause this inflection, 
and is very likely the cause of the large variance in the 3— 
5 /mi continua, as noted by Helou et al. ( 2000|) . Because 
of the characteristics of the L' filter, we conclude that the 
red K — L color must represent a continuum property, 
rather than short- wavelength AFEs. Such a continuum is 
typically accompanied by a lower F6.75/-F15 ratio (see Fig. 
||), due either to AFE suppression or to an F15 excess (or 
both). Either way, we are observing a significant fraction 
of hot dust in the central regions of these galaxies, with 
physical conditions that favor red K — L, namely strong 
Go and high pressure. 



5.1. Physical conditions in star forming regions with 
hot dust 



Malhotra et al. ( 2001 ) describe a physical picture in which 
the infrared line emission arises in a PDR which surrounds 
an expanding H 11 region. They suggest that the ISO diag- 
nostics are probing the typical distance from an OB star 
or cluster to the PDR gas, rather than the global rate of 
star formation. Our new 4 /im data are consistent with this 
picture. Red K — L tends to be associated with an intense 
UV radiation field Go- The trend shown in Fig. || suggests 
that rising 4 /im continua are observed when UV radiation 
field Go ~ 10 3 ' 5 . It is easy to show that such an intense 
Go must be found within a few tens of pc from the ionizing 
star cluster. Assuming the FUV luminosity of a single star 
to be 10 39 erg s _1 (!), and with a star cluster consisting of 
100 OB stars, the typical distance for Go = 10 3,5 is 10 pc. 
In the same conditions, at a distance of 100 pc we would 
need 10000 massive stars, more than all but those in the 
most luminous Super Star Clusters (SSCs: e.g., Calzetti et 
al. |1997| , Turner et al. |2000| ). With the lower FUV lumi- 
nosity from lower-mass stars, the distance to maintain the 
intense Go would be even smaller. Hence, it is likely that 
red K — Lis probing a region within a few tens of pes from 
star clusters containing a few hundred OB stars. Such sizes 
are more than a hundred times smaller than the projected 
size of our observing aperture at the median distance of 
the sample (see Sect. ||). 

Hot dust as measured by high-resolution 10-12 /im im- 
ages of infrared-luminous starbursts also tends to be com- 
pact, and compact nuclear starbursts generally dominate 



their starburst activity (Soifer et al. 2001). Typical MIR 
sizes in the Soifer et al. sample range from < 125 pc to 
~ 1 kpc; however these may be partly overestimated since 
the spatial resolution for the nearest galaxy in their study 
(NGC 3690) is 210pcarcsec _1 which procured mostly up- 
per limits (< 125 pc) for the size of the emitting regions. 
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The hot dust that gives rise to red K — L must be 
heated by massive stars, and according to the PDR mod- 
els, tends to be found in high-pressure environments and, 
from the above discussion, in close proximity of the H II 
region. This intense environment suppresses AFE emis- 
sion, probably because of the strong UV radiation field, 
as proposed by earlier work (e.g., Normand et al. 1995| ). 
In the framework of an expanding H n region surrounded 
by a PDR, such conditions are likely to obtain in young 
H II regions, still embedded within their natal molecular 
cloud, which have not yet had time to break through the 
cloud surface. At first glance then, K — L appears to be a 
measure of age: hot dust is found preferentially in young 
H II regions. 

Nevertheless, the attribution of the rising 4/im con- 
tinuum to youth has some problems. Some very young 
systems have "normal" photospheric+extinction K — L 
colors: the two brightest star clusters in NGC 1569 contain 
Wolf-Rayet (W-R) stars and t heir age has been estimated 
at < 5 Myr (Origlia et al. |200lj) . We measure for NGC 1569 
K — L = 0.7, consistent with stellar photospheres + ex- 
tinction, not with hot dust. Haro2 (not in this sample) 
also contains W-R stars (Vacca & Conti 1992), and as 
such must be younger than 5 Myr, but K — L = 0.5 (Thuan 



19831 ) . Thus, young systems are not necessarily associated 
with red K — L and hot dust. 

It also might be that K — L is a diluted version of 
the flux ratio FIR/B, since variations in FIR/B, or star- 
formation activity, might be translated into variations of 
the hot-dust versus stellar fraction in our observing aper- 
ture. We investigate this in the top panel of Fig. fj], where 
we have plotted Log(FIR/ B) vs. K — L. There is no cor- 
relation between these two quantities, although the galaxy 
with the highest value of FIR/B is associated with the 
reddest K — L (IRAS 23365+3604). It is thus difficult to 
interpret K — L as an indicator of generic star formation, 
since the hot dust traced by K — L appears to be inde- 
pendent of the cooler grains traced by FIR. 



5.2. "Active" and "passive" star formation 

We therefore propose an alternative explanation, namely 
two "extreme" distinct modes of star formation. 
Recent models have shown that dust heating and 
molecular hydrogen production are more efficient in 
dense ( ^ 500 cm" 3 ), compact ( 100 pc) environments 



(Hirashita et al. 2002 ). Moreover, in low-metallicity BCDs, 
the size of the brightest star-forming complex and its ion- 
ized gas density are very well correlated, with the densest 
regions being also the most compact (Hunt et al. 2002| ). 
In this same BCDs sample, size and density are indepen- 
dent of age, since ages derived from recombination line 
equivalent widths arc all £ 5 Myr. The correlations in 
the KP sample studied here seem to suggest a similar 
"dichotomy" since a strong UV radiation field, compact 
size, and pressure are associated with hot dust. The most 
compact dense regions in the BCD study by Hunt et al. 
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Fig. 7. Top panel: Log of the far-infrared-to-blue flux 
ratio FIR/B vs. K — L. FIR/B was taken from Dale et 
al. (2000), with FIR calculated according to the canonical 
depen dence on Fqq and Fioo (Lonsdale Persson & Helou 
1987 ). Bottom panel: Log Li?/^ /area versus K — L. Lfir 



is the far-infrared luminosity and the area is the one of the 
ellipse with axes of the size reported in NED; L fir /area is 
in ergs -1 kpc -2 . Galaxies with corrected H — K < 0.35 
are shown as open circles, and H — K > 0.35 as filled 
ones. The vertical dotted line separates the "active" and 
"passive" regimes (K — L 1). 



(2002) were interpreted by them as indicative of "active" 
star formation, characterized by an intense physical envi- 
ronment harboring several hundreds/thousands of massive 
stars, and capable of efficiently heating dust. At a compa- 
rable star-formation rate, the less dense, less compact H n 
regions, termed "passive" by Hunt et al. (2002), arc not 



able to heat dust efficiently bec ause o f the more diluted 
radiation field (Hirashita et al. 2002 ). We would argue 
that a rising Afim continuum is yet another signature of 
"active" star formation. A falling 4 /im continuum would 
be a sign of "passive" star formation, since the physical 
conditions are not able to produce a measurable fraction 
( ~ 5%) of hot dust; these more diffuse regions tend to 
have lower pressures, a less intense radiation field, and 
prominent AFEs in the MIR spectrum. 

It might be argued that the two extremes are simply 
a measure of the star- formation rate (SFR). We have at- 
tempted to test this hypothesis by plotting the far-infrared 
luminosity per unit area Lpm/a^rea against K — L, since 
Lfir is frequently used to measure SFRs (e.g., Thronson 
& Telesco p.986| ). The bottom panel of Fig. | shows that 
the two quantities are not significantly correlated, nor is 
K — L correlated with Lpm, that is without normalizing 
by the area (not shown). We already said that, although 
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aperture effects may be important because of the large 
IRAS apertures, the top panel of the same figure shows 
no correlation between K — L and the L pjr normalized to 
the B band, nor with F 60 /F W0 (Sect. 3.1), yet another star 



formation indicator of common use. It appears, therefore, 
that the difference between the two extremes is of quali- 
tative rather than quantitative origin. 

It should be emphasized that "active" and "passive" 
star formation as we are defining them are local concepts: 
in a single galaxy, different regions may be character- 
ized by either mode. In our Galaxy AFEs are notably 



diminished in H11 regions (Cesarsky et al. 1998), and the 
MIR continuu m tend s to be stronger in denser structures 
(Abergel et al. 2002 ). M82 is an another example, since 
where the infrared continuum peaks, there are few if any 
AFEs, but at the H 11 region- mol ecular cloud interface, the 
AFEs dominate (Normand et al. 1995 ). On galactic scales, 
there is also evidence that AFEs are diminished as the in- 
tensity of star formation increases (Forster Schreiber et al. 
2002 ). With our ground-based data, we are investigating 
the properties of the dominant mode of star formation in 
the observing aperture, weighted by the brightness of the 
components. At longer wavelengths, higher spatial reso- 
lution is necessary to separate the two modes because in 
large apertures (e.g., IRAS, ISO), the cooler "standard" 
ISM virtually overwhelms any hot dust emission. 

We further note that the bulk of star formation associ- 
ated with starbursts is not necessarily of the "active" type. 
AFEs are always detected in galaxy nuclei with H 11- region 
spectra (Roche et al. 1991); indeed these features domi- 
nate the integrated spectra of prototypical starbursts such 
as M82 and NGC253 (Sturm et al. |2000[ ). They are also 
used to distinguish star formation from nuclear activity as 
the mechanism powering ultra -luminous infrared galaxies 
(ULIRGs, Genzel et al. 1998 ). This apparent "paradox" 
can be understood through the association of AFE s with 
the cool dust phase of the PDR (Haas et al. |2002|) . The 
hot dust we measure with K — L is heated by the intense 
Go close to the Hn region, and without high-resolution 
( ^ 100 pc) measurements, the copious emission from the 
diffuse cooler dust and AFEs in the PDRs dominates the 
observed flux. 

We have found evidence for hot dust in roughly 1 /3 of 
the KP galaxies observed. It is not clear, however, whether 
the properties of these "active" star-forming galaxies are 
"extreme" compared to other examples of nuclear star- 
bursts. An upper limit to the bolometric surface brightness 
of thes e ob jects can be deduced from the estimates given 
in Sec. 



5.1 



if we have a nuclear star cluster containing 100- 
1000 massive stars concentrated in a region of 10-100 pc 
in diameter, the bolometric surface brightness would be 
roughly 10 11 — 10 12 Lokpc -2 . This coincides approxi- 
mately with the maximum mean surface brightness in the 
sample of low- and high-redshift starbursts of (Meurer et 
al. 1997). However, our estimate refers to the nuclear star 
clusters, not to the mean surface brightness over the ob- 
serving aperture as in Meurer et al. As star clusters go, 
those we have detected with K — L are not extraordinary, 



since they are several times less bright than the resolved 
UV clusters described by Meurer et al. ( 1995| ). 

What we have called the "active" mode is character- 
ized by the formation of compact star-forming complexes 
within a high- pressure environment in an intense Go . High 
pressure is also implicated in the strong galactic winds 
common in starbursts (Heckman et al. 1990). Such con- 



ditions also favor the formation of dusty SSCs (Bekki & 
Couch [2001 ) , which are frequently found in merging galax- 
ies (e.g., the Antenn ae: Whitmore et al. [1995|; NG C 1741: 
John son et al. |1999| ), ULIRGs (Scoville et al. |200C| , Shioya 
et al. 2001]) , an d nuclear starbursts (e.g., M82: Gallagher 



& Smith |1999| ; NGC 253: Keto et al. |1999| ; NGC 4214: 
Leitherer et al. 1996| ). Although somewhat uncertain, the 



nuclear star clusters inferred from our data have proper- 
ties similar to SSCs, which could be yet another signature 
of "active" star formation. 



5.3. Speculations for high-redshift studies 

The physical environment of the "active" mode, especially 
high pressure, results naturally from interactions or merg- 
ers (Bekki & Couch 2001). Interestingly, interactions and 

L (Joseph et al. 



mergers also produce red K — L (Joseph et al. 1984). 
According to the Cold Dark Matter hierarchical cluster- 
ing scenario of galaxy formation, the frequency and in- 
tensity of interactions are expected to increase with red- 
shift, and semi-analytic recipes for galaxy formation based 
on "collisional starbursts" are successful at reproducing 
SFRs, colors, mo rphol ogy, and mass trends with redshift 
(Somerville et al. 2001). Presently, the bulk of star forma- 
tion in galaxies, including starburst galaxies and ULIRGs, 
takes place in the "passive" regime. But, due to the higher 
probability of strong tidal stresses and mergers, this situa- 
tion may have been reversed in the past with the "active" 
regime being the preferred way to form stars at high red- 
shift. 

In this context, it is puzzling that we maintain that 
ULIRGs host passive star formation, although most if not 
all of them are mergers. The most plausible explanation 
of this paradox involves spatial resolution and the rela- 
tive dominance of the "passive" ISM in ULIRGs (see Sect. 
|5.2| ). Given the large contribution from their cool ISM, the 
emission in ULIRGs integrated over large spatial scales ap- 
pears passive. Even though "active" star formation must 
be occurring, it can be revealed only on small spatial 
scales, or by its hot dust. A case in point is the ULIRG in 
our sample, IRAS 23365+3604. This galaxy hosts several 
young (<10Myr) star clusters (Surace et al. ^000[) , but 



is classified optically as a LINER (Veilleux et al. [1995). 
The LINER classification is probably not caused by an 
AGN, but rather by sho cks generated in galactic super- 



winds (Lutz et al. 1999 ). We classify this object as "ac- 



tive" because of its extremely red K — but its ISO 



The red K — L cannot easily be attributed to nuclear ac- 
tivity, since LINERs tend to have NIR co lors d ominated by 
normal stellar populations (Lawrence et al. 1985). 
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spectrum contains prominent AFEs (Tran et al. 2001). 
The young star clusters probably heat the dust responsi- 
ble for the red K — L, while the ISO aperture encompasses 
the entire galactic disk (see Surace et al. 2000). Hence, 



IRAS 23365+3604 appears both active and passive: active 
because of its red K — L, and passive because of the over- 
whelming disk ISM. 

Including the effects of active star formation in mod- 
els of galaxy evolution (e.g., Granato et al. 2000) would 
not be straightforward. Nevertheless, when star formation 
occurs on spatial scales lOOpc (see Sect. 5.1), the in- 



tense radiation field, high pressure, and diminished AFEs 
should be taken into account. Indeed, suppressed AFEs 
coupled with the hot dust revealed by our photometry 
substantially alter the infrared spectrum of star-forming 
galaxies. This may be especially important in hierarchi- 
cal merger models because of the supposed link between 
"active" star formation and interactions. 

6. Conclusions 

With our JHKL 1 photometry, we have analyzed the 4 /im 
continuum and its relation with the MIR spectrum and 
FIR emission lines. We find the following: 

1. The majority of the 26 KP galaxies have a falling 4 fim 
continuum, K — L < 1.0, consistent with stellar pho- 
tospheres and moderate dust extinction. 10 of them 
have a flat or rising 4/zm continuum, K — L ^ 1.0, 
consistent with a measurable fraction of 600-1000 K 
hot dust. 

2. K — L is anticorrelated with ISO ratios F6.75/-F15 and 
IRAS ratio F 12 /F 2 5, but only weakly with [Cii]/[Oi], 
and not at all with [C ii]/FIR/ or IRAS ratio F 60 /Fi o- 

3. PDR models for these galaxies show that the hot dust 
measured by red K — L is associated with high pres- 
sures and intense far- ultraviolet radiation fields in com- 
pact ( 100 pc) regions. 

4. These results taken together suggest that star forma- 
tion in these galaxies occurs in two "extreme forms" : 

(a) a relatively rare "active" mode characterized by hot 
dust, suppressed AFEs, high pressure, intense ul- 
traviolet radiation field, and compact size; 

(b) a more common "passive" mode characterized by 
photosphcric K — L colors, with moderate extinc- 
tion, and less extreme physical conditions. 

5. The physical conditions we infer for the star- forming 
regions containing hot dust are similar to those created 
by interactions and mergers. We speculate that such 
intense episodes may have been more common in the 
past, so that the "active" regime could dominate star 
formation at high redshift. 
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